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I We report results of a search for Weakly Interacting Massive Particles (WIMPs) with the sil- 

• • icon detectors of the CDMS II experiment. A blind analysis of 140.2 kg-days of data revealed 

^ three WIMP-candidate events with an expected total background of 0.7 events. The probability 

that the known backgrounds would produce three or more events in the signal region is 5.4%. A 

profile likelihood ratio test of the three events that includes the measured recoil energies gives a 
^ 0.19% probability for the known-background-only hypothesis when tested against the alternative 

WIMP+background hypothesis. The highest likelihood occurs for a WIMP mass of 8.6 GeV/c^ and 

WIMP-nucleon cross section of 1.9x10"^^ cm^. 

PACS numbers: 14.80.Ly, 95.35. +d, 95.30.Cq, 95.30.-k, 85.25.0j, 29.40.Wk 



There is now overwhelming evidence that the bulk of 
the matter in our universe is in some nonluminous, non- 
baryonic form [J. Weakly Interacting Massive Particles 
(WIMPs) [7 form a leading class of candidates for this 
dark matter. Particles of this type would be produced 
thermally in the early universe and are predicted by many 
theoretical extensions to the Standard Model of particle 
physics m O S]. If WIMPs do constitute the dark mat- 
ter in our galaxy, they may be detectable through their 



elastic scattering from nuclei in terrestrial particle de- 
tectors [5 . Numerous experimental groups have sought 
to detect such scattering events using a wide variety of 
technologies [6 . 

The Cryogenic Dark Matter Search (CDMS) collabora- 
tion identifies nuclear recoils (including those that would 
occur in WIMP interactions) using semiconductor detec- 
tors operated at 40 mK. These detectors use simulta- 
neous measurements of ionization and non-equilibrium 
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phonons to identify such events among the far more nu- 
merous background of electron recoils. 

The low atomic mass of Si generally makes it a less 
sensitive target for spin-independent WIMP interactions 
relative to the larger coherent enhancement of the scat- 
tering cross section for heavy nuclei. On the other hand, 
the lower atomic mass of Si is advantageous in searches 
for WIMPs of relatively low mass ('^10 GeV/c^) due to 
more favorable scattering kinematics. New particles at 
such masses are generally disfavored in fits of models to 
precision electroweak data {e.g. [7), but viable models 
in this regime do exist {e.g. [8, 9 ). Renewed interest in 
this mass range has been motivated by results from the 
DAM A/LIBRA [10], CoGeNT [11 , and CRESST [12] ex- 
periments, which can be interpreted as evidence of low- 
mass WIMP scattering. 

During 2003-2008 the collaboration operated CDMS II, 
an array of Ge and Si detectors located at the Soudan Un- 
derground Laboratory [13T[l9]. In its final configuration, 
the CDMS II array consisted of 30 Z-sensitive ionization 
and phonon (ZIP) detectors: 19 Ge (~239 g each) and 11 
Si (~106 g each), for a total of ~4.6 kg of Ge and ~1.2 kg 
of Si. We discriminate nuclear recoils from background 
electron recoils using the ratio of ionization to phonon 
recoil energy (ionization "yield"). Electron recoils that 
occur within ~10 jam of a detector surface can exhibit 
reduced ionization collection. These events are identi- 
fied by phonon pulse-shape discrimination. Our overall 
misidentification rate of electron recoils is less than 1 in 
10^ 

We consider data from the Si detectors using the final 
four run periods of the full CDMS II detector installa- 
tion, acquired between July 2007 and September 2008. 
The Ge results from this data set have been described in 
previous publications [17 . Of the 11 Si detectors, three 
were excluded from the WIMP-search analysis: two due 
to wiring failures that led to incomplete collection of the 
ionization signal and one due to unstable response on 
one of its four phonon channels. Periods of poor perfor- 
mance, as identified by a series of Kolmogorov-Smirnov 
tests, were also excluded from analysis. After all such 
exclusions, the data collected by the 8 Si detectors con- 
sidered in this analysis represent a total exposure of 140.2 
kg-days prior to the application of the WIMP candidate 
selection criteria. 

The responses of these detectors to electron and nu- 
clear recoils were calibrated using events from extensive 
exposures to ^^^Ba and ^^^Cf sources in situ at Soudan. 
Electron recoils from the former were used to empirically 
characterize and correct for the dependence of phonon 
pulse shape on event position and energy. The 356 keV 
gamma ray from the ^^^Ba source has a ~ 4.2 cm atten- 
uation length in Si, and thus the Si detectors generally 
do not show a clear line at 356 keV. Their energy scales 
were calibrated using 356 keV events with total energies 
shared between the Si detector and a neighboring detec- 
tor. 

WIMP-candidate events were identified by a series 
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FIG. 1. Efficiency as function of recoil energy after appli- 
cation of each WIMP-selection criterion shown. Each curve 
from top to bottom shows the cumulative effect of succes- 
sive cuts on the data, such that the bold solid curve shows 
the overall efficiency of this analysis. The abrupt drops in 
acceptance at low recoil energies reflect the elevated energy 
thresholds chosen for some detectors. 



of selection criteria. Events in and near the WIMP- 
candidate region were automatically removed from the 
data set during the analysis, and all WIMP-selection cri- 
teria were defined blindly using calibration and remaining 
WIMP-search data. Thus, WIMP candidates had no im- 
pact on the definition of the selection criteria. A WIMP 
candidate was required to have phonon and ionization 
signals above the noise in exactly one ZIP detector and 
to exhibit no coincident energy in the scintillating veto 
shield. Events in coincidence with the NuMI beam [20 
were also vetoed. We demanded that any candidate event 
occur within the detector's fiducial volume, defined by 
requiring signal consistent with noise in the outer ioniza- 
tion electrode. The recoil energy of each candidate event 
had to lie below 100 keV and above a detector-dependent 
threshold ranging from 7 to 30 keV, chosen blindly based 
on calibration data. Candidate events were further re- 
quired to lie above a detector and dataset-dependent ion- 
ization threshold, which was 4.5<j from the noise as mea- 
sured by randomly acquired triggers. In yield, events 
were required to be within +1.2cr and — 1.8cr from the 
mean of the nuclear recoil yield. Candidate events were 
also required to have phonon pulse timing consistent with 
a nuclear recoil. In order to take advantage of the fact 
that the timing parameters are better measured at high 
energies, the phonon timing data-selection cut was opti- 
mized in three energy bins: 7-20 keV, 20-30 keV, and 
30-100 keV [21]. Fig. [l] shows the "WIMP efficiency," 
i.e., the estimated fraction of WIMP recoils at a given 
energy that would be accepted by these signal criteria. 
The abrupt changes in efficiency are due to the different 
detector thresholds and changes to the timing cuts in the 
three energy bins. Signal acceptance was measured us- 
ing nuclear recoils from ^^^Cf calibration, and is ~40% at 
most recoil energies. After applying all selection criteria. 
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FIG. 2. Ionization yield versus recoil energy in all detectors 
included in this analysis for events passing all signal criteria 
except (top) and including (bottom) the phonon timing crite- 
rion. The curved black lines indicate the signal region {-1.8a 
and +1.2(7 from the mean nuclear recoil yield) between 7 and 
100 keV recoil energies, while the gray band shows the range 
of charge thresholds. Electron recoils in the detector bulk 
have yield near unity. The data are colored to indicate recoil 
energy ranges (dark to hght) of 7-20, 20-30, and 30-100 keV 
to aid the interpretation of Fig. [3] 



the exposure of this analysis is equivalent to 23.4 kg-days 
over a recoil energy range of 7-100 keV for a WIMP of 
mass 10 GeV/c^. 

Neutrons from cosmogenic or radioactive processes 
can produce nuclear recoils that are indistinguishable 
from those from an incident WIMP. Simulations of the 
rates and energy distributions of these processes using 
GEANT4 [22 lead us to expect < 0.13 false candidate 
events (90% confidence level) in the Si detectors from 
neutrons in this exposure. 

A greater source of background is the misidentifica- 
tion of surface electron recoils, which may suffer from re- 
duced ionization yield and thus contribute events to the 
WIMP-candidate region; these events are termed "leak- 
age events". Prior to looking at the WIMP-candidate 
region (unblinding), the expected leakage was estimated 
using the rate of single scatter events with yields con- 
sistent with nuclear recoils from a previously unblinded 
dataset ^\ and the rejection performance of the timing 
cut measured on low-yield multiple-scatter events from 
-^^^Ba calibration data. Two detectors used in this anal- 
ysis were located at the end of detector stacks, so scatters 
on their outer faces could not be tagged as multiple scat- 
ters. The rate of surface events on the outer faces of these 
two detectors were estimated using their single-scatter 
rates from a previously unblinded dataset presented in 
[23] and the multiples-singles ratio on the interior de- 
tectors. The final pre-unblinding estimate for misidenti- 
fied surface electron-recoil event leakage into the signal 
band in the eight Si detectors was OATtoiristat.) events. 
This initial leakage estimate informed the decision to un- 
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FIG. 3. Normalized ionization yield (standard deviations 
from the nuclear recoil band centroid) versus normalized 
phonon timing parameter (normalized such that the median 
of the surface event calibration sample is at -1 and the cut 
position is at 0) for events in all detectors from the WIMP- 
search data set passing all other selection criteria. The black 
box indicates the WIMP candidate selection region. The data 
are colored to indicate recoil energy ranges (dark to light) of 
7-20, 20-30, and 30-100 keV. The thin red curves on the bot- 
tom and right axes are the histograms of the data, while the 
thicker green curves are the histograms of nuclear recoils from 
^^^Cf calibration data. 



blind. After unblinding, we developed a Bayesian es- 
timate of the rate of misidentified surface events based 
upon the performance of the phonon timing cut mea- 
sured using events near the WIMP-search signal region 
[23^ . Multiple-scatter events below the electron-recoil 
ionization-yield region from both ^^^Ba calibration and 
the WIMP-search data were used as inputs to this model. 
Because the WIMP-search sample is sparser compared 
to the calibration data, the combined estimates are more 
heavily weighted towards the calibration data leakage es- 
timates. Additionally the leakage estimate is corrected 
for the fact that the passage fraction of singles and mul- 
tiples differs by a factor of 1.7lo!6' measured on low- 
yield events outside of the nuclear recoil band. The sys- 
tematic uncertainty on the leakage estimate comes from 
the uncertainty on this scale factor, the choice of prior in 
the Bayesian analysis, and the method used to reweigh 
the energy distribution of surface events from calibration 
data to reflect the distribution in WIMP search data. 
The final model predicts an updated surface-event leak- 
age estimate of OAlt.Q'^^{stat.)t.o'^4{syst.) misidentified 
surface electron-recoil events in the eight Si detectors. 
Classical confidence intervals provided similar estimates 
124,. 

After all WIMP-selection criteria were defined, the sig- 
nal regions of the Si detectors were unblinded. Three 
WIMP-candidate events were observed, with recoil en- 
ergies of 8.2, 9.5, and 12.3 keV, on March 14, July 1, 
and September 6 of 2008, respectively. Two events were 
observed in Detector 3 of Tower 4, and the third was ob- 
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served in Detector 3 of Tower 5. These detectors were 
near the middle of their respective tower stacks. Fig. |2] 
ihustrates the distribution of events in and near the sig- 
nal region of the WIMP-search data set before (top) and 
after {bottom) application of the phonon timing criterion. 
Fig. [3] shows an alternate view of these events, expressed 
in "normalized" versions of yield and timing that are 
transformed so that the WIMP acceptance regions of all 
detectors coincide. 

After unblinding, extensive checks of the three candi- 
date events revealed no data quality or analysis issues 
that would invalidate them as WIMP candidates. The 
signal-to-noise on the ionization channel for the three 
events (ordered in increasing recoil energy) was measured 
to be 6.7(7, 4.9(7, and 5.1(7. A study on possible leakage 
into the signal band due to ^^^Pb recoils from ^^^Po de- 
cays found the expected leakage to be negligible with 
an upper limit of < 0.08 events at the 90% confidence 
level. The energy distribution of the ^^^Pb background 
was constructed using events in which a coincident a par- 
ticle was detected in a detector adjacent to one of the 8 
Si detectors used in this analysis. 

This result constrains the available parameter space 
of WIMP dark matter models. We compute upper lim- 
its on the WIMP-nucleon scattering cross section using 
Yellin's optimum interval method [25]. We assume a 
WIMP mass density of 0.3 GeV/c^/cm^, a most probable 
WIMP velocity with respect to the galaxy of 220 km/s, 
a mean circular velocity of Earth with respect to the 
galactic center of 232 km/s, a galactic escape velocity of 
544 km/s [26 , and the Helm form factor [27]. Fig. [I] 
shows the derived upper limits on the spin-independent 
WIMP-nucleon scattering cross section at the 90% con- 
fidence level (C.L.) from this analysis and a selection of 
other recent results. The present data set an upper limit 
of 2.4 X 10-^^ cm^ for a WIMP of mass 10 GeV/c^. We 
are completing the calibration of the nuclear recoil energy 
scale using the Si-neutron elastic scattering resonant fea- 
ture in the ^^^Cf exposures. This study indicates that our 
reconstructed energy may be 10% lower than the true re- 
coil energy, which would weaken the upper limit slightly. 
Below 20 GeV/c^ the change is well approximated by 
shifting the limits parallel to the mass axis by ~ 7%. In 
addition, neutron calibration multiple scattering effects 
improve the response to WIMPs by shifting the upper 
limit down parallel to the cross-section axis by ~ 5%. 

A model of our known backgrounds, including both 
energy and expected rate distributions, was constructed 
for each detector and experimental run for each of the 
three backgrounds considered: surface electron recoils, 
neutron backgrounds, and ^^^Pb recoils. Simulations of 
our background model yield a 5.4% probability of a sta- 
tistical fiuctuation producing three or more events in our 
signal region. 

This model of our known backgrounds was used to in- 
vestigate the data in the context of a WIMP+background 
hypothesis. We performed a profile likelihood analysis, 
including the event energies, in which the background 
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FIG. 4. Experimental upper limits (90% confidence level) for 
the WIMP-nucleon spin-independent cross section as a func- 
tion of WIMP mass. We show the limit obtained from the 
exposure analyzed in this work alone (blue dotted line), and 
combined with the CDMS II Si data set reported in p3l [28] 
(blue solid line). Also shown are limits from the CDMS 
II Ge standard 17 and low-threshold [29 analysis (dark 
and light dashed red), EDELWEISS low-threshold [30^ (long- 
dashed orange), XENONIO S2-only [3T| (dash-dotted green), 
and XENONIOO 32 (long-dash-dotted' green). The filled re- 
gions identify possible signal regions associated with data 
from CoGeNT 33 (dashed yellow, 90% G.L.), DAMA/LIBRA 
[T0l[34i (dotted tan, 99.7% G.L.), and GRESST [12, 35^ (dash- 
dotted pink, 95.45% G.L.) experiments. 68% and 90% G.L. 
contours for a possible signal from these data are shown in 
light blue. The blue dot shows the maximum likelihood point 
at (8.6 GeV/c^ 1.9 x 10"^^ cm^). 



rates were treated as nuisance parameters and the WIMP 
mass and cross section were the parameters of interest. 
We profiled over probability distribution functions of the 
rate for each of our known backgrounds. The highest like- 
lihood was found for a WIMP mass of 8.6 GeV/c^ and 
a WIMP-nucleon cross section of 1.9x10"^^ cm^. The 
goodness-of-fit test of this WIMP+background hypoth- 
esis results in a p- value of 68%, while the background- 
only hypothesis fits the data with a p- value of 4.5%. 
A profile likelihood ratio test finds that the data favor 
the WIMP+background hypothesis over our background- 
only hypothesis with a p- value of 0.19%. Though this 
result favors a WIMP interpretation over the known- 
background-only hypothesis, we do not believe this result 
rises to the level of a discovery. 

Fig. |4] shows the resulting best-fit region from this 
analysis (68% and 90% confidence level contours) on 
the WIMP-nucleon cross-section vs. WIMP mass plane. 
The 90% C.L. exclusion regions from CDMS IPs Ge 
and Si analyses and EDELWEISS low-threshold analy- 
sis cover part of this best-fit region, but the results are 
overall statistically compatible. There is much stronger 
tension with the upper limits from the XENONIO and 



5 



XENONIOO experiments, which significantly constrain 
this parameter space under standard assumptions about 
the WIMP velocity distribution and WIMP-nucleus in- 
teractions. 
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